Semiclassical degeneracies and ordering for highly frustrated magnets in a field 
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We discuss ground state selection by quantum fluctuations in frustrated magnets in a strong 
magnetic field. We show that there exist dynamical symmetries - one a generalisation of Henley's 
gauge-like symmetry for collinear spins, the other the quantum relict of non-collinear weathervane 
modes - which ensure a partial survival of the classical degeneracies. We illustrate these for the case 
of the kagome magnet, where we find zero-point energy differences to be rather small everywhere 
except near the collinear 'up-up-down' configurations, where there is rotational but not translational 
symmetry breaking. In the effective Hamiltonian, we demonstrate the presence of a term sensitive 
to a topological 'flux'. We discuss the connection of such problems to gauge theories by casting the 
frustrated lattices as medial lattices of appropriately chosen simplex lattices, and in particular we 
show how the magnetic field can be used to tune the physical sector of the resulting gauge theories. 



I. INTRODUCTION 

The rich behaviour presented by geometrically frus- 
trated magnets has made them a popular subject of 
study. 1 Particular attention has been devoted to systems 
in which the spins reside on corner-sharing units: the 
kagome lattice in d = 2 dimensions, consisting of corner- 
sharing triangles, and the d = 3 pyrochlore lattice, made 
up of corner-sharing tetrahedra, whose sparse connectiv- 
ity makes them particularly strongly frustrated^ 

Classical ground state degeneracies can be considered 
to be the hallmark of such geometrical frustration. One 
natural question is what happens to these ground state 
degeneracies in the presence of quantum fluctuations. 
Generically, one would expect the degeneracies to be 
lifted, so that they would show up in a quantum model as 
a large density of states at low energies, on a scale set by 
the strength of the quantum fluctuations. Indeed, such 
a large density of low-energy excitations is perhaps the 
most striking feature of the kagome spin S = 1/2 Heisen- 
berg magnet^ although as yet there exists no agreement 
on its origin. 

Given the difficulties in studying the limit of strong 
quantum fluctuations, small S, it has proven to be fruit- 
ful to consider instead the semiclassical limit, 3 S — > oo, 
where quantum fluctuations endow the classical ground 
states with an effective energy functional via the zero- 
point energy of the spin excitations, %0({S^ })/2, 
where the <I> denote the frequencies of the excitations 
around a ground state configuration {S^}. (In the 
following, we set Ti = 1.) Perhaps the most striking 
result was that a sizeable, albeit reduced, degeneracy 
can survive in the semiclassical limit. The most promi- 
nent example of this was again provided by the kagome 
lattice, where the equations of motion for excitations 
around any coplanar ground state can be written in an 
identical form, so that all these states have the same 
zero-point energy^ The next important step was Hen- 
ley's observation; 5 sketched below, that degeneracies for 
collinear ground states could be understood in terms of 



a gauge-like symmetry in the equations of motion. 

In a related development, a great deal of attention has 
recently been focused on the behaviour of frustrated mag- 
nets in an external field, in particular on magnetisation 
plateaux and magnetocaloricsiSiLSiSiifliiiii^iH 

In this paper, we follow up on those two strands of 
work by considering the properties of frustrated magnets 
on corner-sharing geometries in an external field in the 
semiclassical limit. We start from the observation that 
such corner-sharing lattices can be thought of as medial 
lattices of a 'simplex lattice'; the spins are bond variables 
on the simplex lattice, whereas the corner-sharing units 
reside on its sites. 

The ground state condition imposed by classical 
Hamiltonian then takes the form of a local constraint, 
which can conveniently encoded in a gauge theory, the 
'physical' sector of which can be changed by the appli- 
cation of a magnetic field, so that one can tune between, 
e.g., loop and dimer models in this fashion. 

We then consider the harmonic theory of excitations 
in such magnets by generalising the equations of motion 
derived in Ref. |2| to include a magnetic field. This will 
allow us to show that Henley's gauge-like symmetry for 
collinear spins can be generalised to this case, but in 
addition that a different dynamical symmetry exists for 
non-collinear spin configurations. 

We illustrate these ideas with the example of a 
Heisenberg magnet on the kagome lattice in an ex- 
ternal field. This choice allows us to build on a 
substantial body of previous work (on classical equa- 
tions of motion, order by disorder and magnetisation 
plateauxj2ii2iiiii2ii£iiiiiii).In keeping with received wis- 
dom, the semiclassical order by disorder we find favours 
coplanar or collinear states. In different field ranges, both 
states derived from Q = and Q = -s/3 x y/3 configura- 
tions occur; we find Z$ and Z% x U(l) order parame- 
ters, some of which will only be expressed at minuscule 
temperatures. For the value of the magnetic field which 
permits a collinear configuration, we find that the semi- 
classical ordering pattern differs from that proposed for 
S = 1/2M. 
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We contrast this behaviour to the kagome lattice 
thought of as a network of corner-sharing hexagons, 
where the effective semiclassical model does not permit 
such simple solutions. Finally, we comment on experi- 
ment, and close with some more general remarks about 
gauge theories for this problem. 

The remainder of this paper is organised as follows. 
In Sect, m we introduce notation by writing down the 
Hamiltonian in a field for a simplex lattice and deriving 
fSect. mUl the equations of motion. Sects . HVI and IVI deal 
with the dynamical symmetries and the properties of ef- 
fective Hamiltonians. Sect. IVll explores in detail ground 
state selection in the kagome magnet, and addresses is- 
sues arising for the kagome-hexagon model. We conclude 
inSect.rvttl 





FIG. 1: (colour online) The kagome lattice (fat lines) is the 
medial lattice of both the honeycomb and the triangular lat- 
tices (thin lines). The former corresponds to a lattice of 
corner-sharing triangles, while the latter consists of corner- 
sharing hexagons. 



II. MEDIAL LATTICE CONSTRUCTION 

A frustrated lattice made up of corner-sharing units 
(also known as simplices), can be obtained by considering 
any lattice A, which we will refer to as the simplex lattice. 
The midpoints of the bonds of the simplex lattice define 
its medial lattice, T. 

The choices for A of the square, honeycomb and dia- 
mond lattices yield the well-known checkerboard, kagome 
and pyrochlore lattices as their respective medial coun- 
terparts T. In fact, there exists a second simplex lattice, 
the triangular lattice, which has the kagome lattice as its 
medial lattice T (see Fig.^). As discussed further down, 
systems on medial lattices generally lend themselves to 
a ga uge theore tic description, for concrete instances see 
Refs. ll8ll9l2nl 

Since sites i e T are also links a/3 (with a, (3 € A) on 
the simplex lattice, we can alternatively denote the spin 
Si by a link variable, S a p. A brief remark on notation: 
we denote the location of a spin as subscripts, its compo- 
nents as superscripts, and the spin length by S (without 
super- or subscripts). 

The ground state condition of the frustrated spin 
Hamiltonian now imposes constraints on the values of the 
link variables S a p emanating from a site a S A. For in- 
stance, denoting the set of links (i.e., sites of T) by v(a), 
the nearest neighbour exchange Hamiltonian is given by 



H, 



do) 



J 



(2.1) 



Si 



where L a = 

The ground state condition L a = Vet represents such 
a local constraint. An external magnetic field provides a 
handle for tuning this constraint. Adding a Zeeman term 
(where g is the g-factor, /is the Bohr magneton and B 
the magnetic field), 



H z = 



gH B Si ■ B 



(2.2) 



yields the Hamiltonian, 

H = H J + H z = ^(Ll- gj ^B.L a) 

a 

J J^ 2 



JS 2 \ " /2 
Z^'<* 



(2.3) 



where h = gjisB / J , and where we have dropped an (h- 
dependent) constant. All states with l a = L a — h/2 = 
Va are thus ground states, provided such states exist. 
Here and in the following, expressions for spins, fields 
and frequencies are rescaled by the spin length S as in 
L a = L a / S . 

For simplicity, let us first consider the case of Ising 
spins (n = 1, S = ±1) for which the L a are integers be- 
tween ±z a , where z a is the coordination of site a G A. 
Effectively collinear Ising spin configurations can also 
arise in highly frustrated continuous spin models, via dif- 
ferent forms of order by disorder ^LifiiSiiSSiSi If z [ s even 
(odd), so are the possible L a . 

If we denote spins pointing up by a link occupied by 
a dimer, and spins pointing down by an empty link, one 
observes that L = z corresponds to the non-degenerate 
saturated state with all spins aligned, and all links empty. 
As has been noted in many contexts, L = z — 2 is equiva- 
lent to one link emanating from each site being occupied. 
A magnet in a field on the simplex lattice T is thus equiv- 
alent to the hardcore dimer models on its simplex lattice 
A. Similarly, L = z— 4 corresponds to loop models, which 
are fully packed in the sense that each site is part of a 
loop. A familiar example of this is the pyrochlore magnet 
in zero field. 

For sufficiently large z, this continues for further L > 0. 
We have illustrated these cases in Fig. [3 where L = 
corresponds to the easy-axis kagome model of Ref. I24I 
which from the current perspective is a three dimer model 
on the triangular lattice. 

For Ising spins, at the transition field where spin L = 
z—2n becomes degenerate with L = z— 2n+2, the allowed 
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L=4 



FIG. 2: (colour online) Multi-dimer configuration correspond- 
ing to different values of the hexagon magnetisation, L. Be- 
sides the usual dimer model below saturation (L = 4), the 
L — 2 sector defines a loop model, and L = a three-dimer 
model. 



vertex configurations are those of both cases. This is in 
particular the case for odd z in zero field, and it goes 
along with a spike in the zero temperature entropy at 
the field strength in question. 



III. EQUATIONS OF MOTION 



IV. DYNAMICAL SYMMETRIES 

A. Collinear ground states 

For collinear spins in zero field, we have L = (L 1 , L 2 , 0) 
and Si — (0, 0, S 3 ). It was pointed out by Henley that the 
form of the equations of motion remains unchanged when 
making the Ising (£ = ±1) gauge-like variable transfor- 
mation 



together with 



Ca L 



Saf3 — * CaSapCfS 



(4.1) 



(4.2) 



for any set of {£}. 

Thus, provided the new state, with the S a p changed in 
sign, is again a ground state of the Hamiltonian i|2.3[l . the 
spectrum of its excitation energies is going to be identical 
to that of the starting state. 

It can be seen that this continues to be the case for the 
equations of motion for the {/} in the presence of a field, 
which can be seen by writing them out explicitly: 



The simple form of the Hamiltonian in terms of the L 
suggests that these variables should be useful quantities 
to use when constructing effective theories. Indeed, the 
thermodynamics of the pyrochlorc Hciscnbcrg antiferro- 
magnet is well captured by formulating a theory in terms 
of the L of a single tetrahedron^ 

Similarly, the L are useful quantities when discussing 
the spectrum of excitations around a classical ground 
state. In zero field, changes of the spin orientation that 
keep L = cost zero energy. Hence, all finite frequency 
modes involve L ^ 0. Indeed, the zero field equations of 
motion for Heisenberg spins around a ground state con- 
figuration {S'^ )} is£ 



La — J S 



sr c(°) 



0:a 



J af3 



X Lf 



(3.1) 



(CJl)-^(CX) 



(4.3) 



It = 



where we take the field along the 3-direction. 

For a non-collinear state, this transformation does not 
work as the bottom member of Eq. 14.31 does not hold in 
that case. This implies that all spins need to point along 
the field axis, as in the ground state L a = h/2 implies 
that not all S 3 <E v(a) can be zero. 

Further down, we show how the kagome magnet at 
1/3 of its saturation field provides an example of such a 
symmetry, including a subextensive degeneracy as well 
as a resulting Z% order parameter. 



Here and in the following, the symbol x denotes a vector 
product, and the sum Ylp- a runs over au Spa € v(a). 
In the following, we drop the superscript as, in our 
harmonic theory, the spins S a /3 will always refer to a 
ground state configuration. 

By the same token, the simple form of Ea. l2.3l suggests 
formulating the equations of motion in terms of 



h/2 



(3.2) 



The resulting form of the equations of motion have a 
similar structure (see Apd. Ia"|| : 



(3.3) 



fi-.a. 



B. Non-collinear ground states 

There exists a different dynamical symmetry for the 
case of non-collinear spins. This can be seen by writing 
the equations of motion in component form: 



f + -I 2 

i 2 - -i 1 

(fill) 



c-3 ;2 



(4.4) 



Ea. 13.31 corresponds to 77 = 1. 

For 77 = —1, these equations retain their form under 
the transformation 



r 



t)ll ; 5, 



1,2 

a/3 



ol,2 



(4.5) 
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Notice that this transformation is distinct from the 
collincar gauge transformation: for collinear states, Z 3 = 
S 1 ' 2 = 0, and the transformation does not generate a new 
spin state. 

This transformation needs to be applied consistently: 
from the bottom member of Eq. 14.41 it follows that a 
change (|4.5|) involving the plaquette /3 requires a con- 
comitant change on those neighbouring a plaquettes for 
which IS^I 7^ 1. If all IS* 3 ^ ^ 1, one therefore ends up 
transforming all the spins in the system, which amounts 
to a global symmetry operation. 

However, if one can find a set of spins the neighbours 
of which all point in the 3-direction, the transformation 
is non-trivial. This is the same condition as that for 
the weathervane modes of the classical kagome magnet 
in zero field^ The local dynamical Z2 symmetry in the 
semiclassical problem is therefore the vestige of the 'acci- 
dental' local £7(1) symmetry present in the classical case, 
and it persists to the case of nonzero fields. 

Here, we have demonstrated that it holds not only for 
the zero-field case, in which the angle between any two 
neighbouring spins is uniform. Below we provide an ex- 
ample of this for the kagome magnet in a field, and dis- 
cuss the order parameter resulting from the ensemble of 
semiclassically degenerate states. 

V. EFFECTIVE HAMILTONIANS 

The basic idea of effective Hamiltonians is to discard as 
much microscopic information as possible, essentially by 
integrating out fluctuations, in the hope that this proce- 
dure will give rise to an expression which is simple (and, 
ideally, local) £ 




FIG. 3: (colour online) Parametrisation of the Q = ground 
states. The fat arrow denotes the magnetic field; its length is 
h/2. cj> denotes the rotation angle of the top pair of spins out 
of the plane. 
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FIG. 4: Zero-point energy per spin, in units of JS, at h/2 — 1 
as a function of 6 (Fig.0. 



A. Selection of collinear configurations 

It has been suggested, first in the seminal work by 
Shender on quantum order by disorder, that quantum 
fluctuations favour collinear spin configurationsi^iiS 3 . in- 
tegrating them out in a semiclassical framework leads, 
crudely, to an effective biquadratic Hamiltonian of the 
form 

Hi = J A S^2(Si-Sj) 2 . (5.1) 

(ij) 

We have tested this assertion in the present context 
by considering the Q — configurations at a magnetic 
field where a collinear state is possible: at h = 2, the 
up-up-down state is possible, in which the spins on two 
sublattices of the kagome lattice point up, and down on 
the third (Fig. [J. 

At h = 2, there is in fact a two-parameter family of 
possible Q — states, as depicted in Fig. In Fig. 0] 
we have evaluated the zero-point energy as a function of 
9, that is to say for coplanar spin configurations. Here, 
9 = corresponds to the collinear state. 



From this picture, two conclusions can be drawn. 
Firstly, the biquadratic form broadly has the correct 
shape, with the minimum and maximum in the right 
places, and a monotonic behaviour in between. Secondly, 
the shape of the curve is in quantitative disagreement 
with the simple biquadratic term. Therefore, when try- 
ing to extract the effective J4, one should proceed in two 
different ways, depending on what information one is in- 
terested in. 

If the question one is asking concerns the overall en- 
ergy scale at which quantum fluctuations start to assert 
themselves, one is basically interested in the difference 
between the highest and the lowest point on that curve. 
By contrast, if one starts from the lowest temperatures 
and asks about the energy cost of small deviations from 
collinearity, one needs to fit the curvature at low 9. (In 
fact, fitting to a power law gives an exponent deviating 
substantially from 2 already for small 9; however, this 
exponent does seem to extrapolate to 2 within the limits 
of our numerical accuracy) . 

Note that the values for J4 thus obtained differ by an 
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order of magnitude when the low-0 fit is done to be ac- 
curate around 9 = 5° (inset of Fig. ■ This discrepancy 
continues to grow as the range of 9 is lowered. 

Is the form of the effective biquadratic interaction more 
accurate in the absence of a field? For the kagome mag- 
net, there are no collinear states available in zero field. 
However, for the checkerboard (planar pyrochlore) mag- 
net, such configurations do exist, and we have checked 
that the corresponding Eq(9) has the same behaviour in 
that case. 



B. Gauge form of Hamiltonian 

By defining Z ± = l x ±il 2 , one can rewrite the equations 
of motion for collinear spins as 

±ult = ^5 a e + S a p)l± . (5.2) 

To obtain Eq = JS(\ui\)/2, this form allows us to make 
use of Henley's trick for calculating the zero point energy 
by writing 

Tr|w| = TrVu^ (5.3) 



= Try ^— + zj 5 a p + hS a/3 + M ap , 

and then Taylor expanding the square root to evaluate 
the trace term by term. Here z is the coordination of 
the simplex lattice, and M a p = S a7 S^p — z5 a p, so that 
its diagonal entries vanish. Thus, formally (and up to a 
constant), for N a p = — (hS a p + M a p) / V 'h 2 / '4 + z, 



Tr|w| = y — + z Tr y/S a p - N a p 

= -\/?^ t (2 rf ™ m (^) 

V 4 ^ 2 m m! v ' 

m=2 

TriV m , as in zero field, only contains terms of the type 
S a pS[}~f...S ua , which can be identified with closed loops 
a — *■ 7 — ► • • • — >£j— > a on the lattice A as each site 
of A occurs an even number of times. 

Besides altering the coefficients of different loops, the 
inclusion of the magnetic field has changed the structure 
of this expression by including a term linear in S a p in the 
Taylor expansion, thus allowing for the possibility that 
TriV™ contains terms with an odd number of factors of 
S a p. For bipartite lattices A, such loops are not possi- 
ble, but they do exist for nonbipartite ones, where the 
structure of the expansion for h =/= Q thus contains terms 
which are absent in zero field. 

In any case, these loops have the form of fluxes $ - the 
factors 

$ = S a pSp~t...S U a (5.5) 

have the appearance of line integrals of the Ising 'gauge 
field' S. Knowing the form of the terms appearing in 



the Hamiltonian, one still needs to fix their coefficients. 
This problem has recently been adresscd in some detail 
by Hizi and Henley^ 

Here, we make some further remarks on the structure 
of the relevant solutions for the case that the dominant 
terms in the effective Hamiltonian are those involving 
the elementary plaquettes of the simplex lattice A. By 
elementary plaquettes we mean the plaquettes the tiling 
of which makes up A, for instance a square in the case of 
the square lattice. 

For a two-dimensional lattice A which is a tiling of 
such elementary plaquettes, fixing the flux through the 
elementary plaquettes determines the fluxes on all con- 
tractible loops of the system. If the pattern of elementary 
fluxes is the same in two different states, this means that 
the contribution to the energy is the same in both states 
even for longer loops. 

This does not mean that the energy is uniquely de- 
termined by specifying the elementary fluxes: the fluxes 
through non-contractible loops such as a loop winding 
around a cylinder in the case of periodic boundary con- 
ditions are not thus determined. Below, we present a nu- 
merical determination, for the case of the kagome lattice, 
that two states having locally the same flux but differing 
in the flux through a non-contractible loop have an zero- 
point energy difference which approaches a constant, JS, 
in the thermodynamic limit. 

Therefore, even though the energy measured locally 
in the flux expansion is identical everywhere, there is a 
nonvanishing nonlocal contribution to the energy. 

VI. GROUND STATE SELECTION IN THE 
KAGOME MAGNET 

In this section, we give examples of the dynamical sym- 
metries discussed above, and discuss the effective models 
and order parameters which arise. 



A. Collinear states 

We start with the case of collinear spins at a magnetic 
field h/2 = 1, where L = 1 allows for states in which two 
spins in each triangle are aligned with the field, and one is 
antialigned. As mentioned above, these states map onto 
dimer model on honeycomb lattice, by placing a dimcr 
on links where S a p points down (Fig.^J). 

There are exponentially many such dimer coverings, 
leading to an entropy of S = 0.1077 per spin once the 
collinear configurations have been selected. The zero- 
point fluctuations will lift the remaining degeneracy as 
different collinear configurations have different 'fluxes' 
through their elementary plaquettes. 

Two particularly simple representatives are the Q = 
configuration (all dimers residing on vertical links, 
Fig. [5}, and the Q = y/3 x y3 configuration, which in- 
volves a tripling of the unit cell. The former has uni- 
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FIG. 5: (colour online) The kagome magnet in a field h = 2. FIG. 6: Energy cost of a topological 'flux' in the collinear 

Down spins are denoted by dimers. A gauge-like transforma- kagome ground state depcited in Fig. as a function of sys- 

tion is carried out on a string of triangles: the spins at the tem size, 
corners of each shaded triangle are flipped. 



formly zero flux in each plaquette, whereas the latter has 
a flux of 7r in two plaquettes out of three. 

We have evaluated the zero-point energy for these 
states for a large system (and for all collinear configu- 
rations of a finite system, see below); as we disagree in 
part with several of the published results on this, we have 
given details of the calculation in App. [Bj The fact that 
the Q = state has a lower zero would imply that fluxless 
configurations are preferred. 

It turns out to be easy to characterise these config- 
urations completely. Only the states in which exactly 
two dimers occur on each hexagon are allowed, for peri- 
odic boundary conditions. For the readers familiar with 
height models, these are the maximally tilted states, or 
equivalently the non-flippable ones. 

For a system of linear length I, with 3/ 2 spins and peri- 
odic boundary conditions, we find the number iVo of zero 
flux states to be 

JV =3(2 , -1). (6.1) 

However, not all of these can be related to one another 
via the transformations described above. 

To see this, let us consider in more detail what kind of 
transformations are generated by 1)4.11 14.2(l . Note that, 
applied to a single triangle a, the total spins of triangle a 
and the neighbouring f3's change, so that these violate the 
ground-state condition. To undo this change, one needs 
to apply the transformation to one-dimensional sets of 
triangles, as depicted in Fig. [SJ This tranformation is 
analogous to the one discussed by Tchernyshyov et al. 
for the planar pyrochlore lattice; the zero-field case dis- 
cussed there exhibits the additional feature of indepen- 
dent transformations on different sublattices^ 

The dimers of the resulting configuration, superim- 
posed on the initial configuration, form two 'strings' 
(lines of alternating new and old dimers) . Even though a 



single string does not change the local fluxes, such strings 
can only be created in pairs by the gauge transformation 
(and they can be separated from one another by further 
transformations) . 

One would therefore in general expect the configura- 
tions differing by an odd number of strings to differ in 
energy; we have checked that this is indeed the case, see 
Fig. the energy cost for an odd number of strings re- 
markably equals, to within our numerical accuracy, pre- 
cisely JS, which suggests that this is a robust result. As 
the fluxes through any contractible loop are the same in- 
dependent of the parity of the number of strings, this 
energy difference in the effective flux Hamiltonian is due 
to a topological flux. 

To find the number of ground states, we thus need 
to enumerate the number of even, iVg, and odd, Nq 
configurations separately. For odd I, Eg. 16.11 gives the 
correct number of ground states as the flux around a 
non-contractible loop Y[ Sap always equals —1 in one 
direction. For / even, we find 

Nl = 3(2' _1 - 1) , N% = 3(2 i " 1 ) . (6.2) 

We have checked these results against numerics: we 
have computed the zero-point energy for all collinear con- 
figurations for / = 4, 5 and 6, the total number of which 
is 417, 7623 and 263640, respectively. The number of 
ground states is given by 21, 93 and 93, respectively; 
these include the Q = states. For the even case, the 
odd first excited states number 24 and 96, in accordance 
with the above expressions. 

What is the order parameter describing this ensemble 
of configurations? In the presence of a field, the SU(2) 
symmetry is broken down to U(l), but this symmetry is 
unbroken for the collinear states. The ensemble, how- 
ever, does break (real space) rotational symmetry. This 
happens because the strings defined above cannot cross, 
so that any configuration only includes dimers of at most 
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two different orientations. Also, shifting the dimers along 
one string (e.g. by having a defect tunnel along the 
string) connects only dimer configurations with the same 
dimer orientation absent. The three pure Q = (zero- 
string) configurations are special in this respect as they 
are the only ones to contain only one dimer orientation. 

The phase at the collinear point therefore is described 
by an order parameter with Z3 symmetry. Note that for 
5=1/2, aQ = \/3 x \/3 ordering has been proposed r 3 - 
so that there should be a phase transition as 5 is lowered. 
This appears reasonable given the fact that at small 5, 
collective quantum tunnelling around hexagons can be 
expected to play a role. These processes, which are ab- 
sent in the 1/5 expansion, will favour the Q = V3 x \/3 
structures i2i 



B. Coplanar states 

The analysis in this case is quite analogous to the pro- 
ceeding subsection. Briefly, we have checked explicitly 
(in intervales of h of 0.2) that, among all Q = con- 
figurations, the favoured ones are the same as found in 
previous studies for triangle-based lattices; these are de- 
picted in Fig. El 7,28 The spin orientations for < h < 2 
are given as 

51 = (0,0,-1) 

5 2 = (sin 6,0, cos 6) (6.3) 

5 3 = (-sin 0,0, cos 6) 
2cos6> = h/2 + 1. 

For 2 < h < 6, one has 

5i = (sin/3, 0, cos/3) 
52 = 53 = (sin a, 0, cos a) (6.4) 
2 cos a = (3 + h 2 /A)/h , 2cos/3 = h - 4 cos a. 

Fig.0shows the zero-point energy of the Q = config- 
uration as a function of h. As for the collinear case, we 
compare the energies of Q = and Q = \/3 x v3; this 
is shown in Fig. |SJ The data in that plot has been ob- 
tained by numerically integrating the zero point energy 
using mathematica, and by extrapolation of the finite size 
results. 

We find that for h < h p , the Q — configuration is 
always energetically favoured over the Q = \/3 x y/3. For 
a region of finite width for h > h p , the Q = V3 x y3 is in 
fact the one with a lower energy. In any case, note that 
the energy differences are very small - the Q = y/3 x 
a/3 is never lower in energy by more than a few 10~ 3 JS 
per spin. The location where the energy difference is 
appreciable is around the collinear point, h = 2. 

Let us now look at the order parameter involved which 
exists in the ensemble of states connected by the gauge 
transformations. 

The vector normal to the plane of the spins defines an 
azimuthal angle with respect to the field direction; this 
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FIG. 7: Zero-point energy per spin of coplanar kagome Q = 
configurations as a function of the magnetic field. 

reflects the residual U(l) symmetry. This symmetry will 
be broken in the coplanar states. 

For h > 2, all spins have 5 3 ^ and hence our dynam- 
ical symmetry (|4.3|l cannot be used to obtain a degener- 
ate configuration. [There do exist degeneracies between 
states at higher energies.] As two of the three spin di- 
rections are the same, this will again mean that only the 
position of one sublattice is special, so that the order 
parameter will be Z3xU(l), on top of any translational 
symmetry breaking. 

By contrast, for ft, < 2, all the three spin directions are 
different. However, the spins on one sublattice points 
along the field direction, and hence applying the trans- 
formation (|4.3|) leads to other degenerate configuarations. 
The strings involved in the dynamical symmetry inter- 
change the spins with 5 3 > and alternating ±5 , so 
that in the overall ensemble, S1.2 have no preferred sub- 
lattice, whereas 53 does. The resulting symmetry break- 
ing is therefore restricted to the sublattice on which the 
spin with 5 3 = — 1 resides. The order parameter is again 
Z 3 xU(l). 



C. Kagome-hexagon model 

In the collinear examples so far, the effective gauge 
Hamiltonian has been simple in that there have existed 
configurations which minimise all the terms for the short- 
est loop that arises. These effective classical Hamiltoni- 
ans are therefore unfrustrated. Different scenarios, how- 
ever, are possible. For instance, in the pyrochlore slice, 29 
the actual effective Hamiltonian is not dominated by the 
shortest loop. 

Here, we display a third way the effective Hamiltonian 
can turn out. This occurs in kagome-hexagon model, in 
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FIG. 8: (colour online) Zero-point energy difference per spin, 
in units of JS, between Q = v3x v3 and Q — obtained from 
a numerical evaluation of the Brioullin zone integral (line) of 
the frequencies, and from an extrapolation of the zero-point 
energies of finite systems (circles). 



which the kagome lattice arises from a triangular sim- 
plex lattice, that is to say, it is considered as a lattice of 
corner-sharing hexagons. The zero-field Ising case has 
been studied in Ref. 0- I n the present situation, as 
discussed in Sect. [V] the collinear manifold is selected 
dynamically, and the first sign of the presence of quan- 
tum fluctuations should be incipient nematic correlations 
(never long-ranged, by Mermin- Wagner) as the temper- 
ature is lowered. 

For the collinear states with L = 4, which can be 
mapped onto the dimer coverings on the triangular lat- 
tice (Fig. EJl, it is the external field which determines the 
direction of the collinear axis, and no spontaneous sym- 
metry breaking is required. 

The effective flux Hamiltonian for the resulting dimer 
model differs from the previous cases in several ways. 
Firstly, loops which appear in the expansion (|5.4I) can 
have odd lengths, as the simplex lattice is non-bipartite. 
In fact, the shortest loop has length three. 

Secondly, the term in the effective Hamiltonian corre- 
sponding to the shortest loop does not act to lift the de- 
generacy between the dimer states: the number of loops 
with an odd number of dimers equals twice the number 
of dimers; all the other loops contain no dimers. This 
follows from the collinear ground-state constraint (the 
hardcore dimer condition) and is independent of dimer 
configuration. 

Thirdly, whereas this first term is trivial, the next 
term, corresponding to loops of length four, is frustrated: 
there is no state in which all loops of length four contain 
an even (or all an odd) number of dimers. The resulting 
effective Hamiltonian is therefore simple in that it is a 
purely classical one, but still complicated, as it is frus- 
trated. We leave an analysis of this model as a subject 
for future work. 

In closing this section, we note that, unlike in the con- 



ventional kagome model, the selection of collinear states 
does not take place classically, by the counting argument 
of Ref. 2, Rather, this is an example where there will be 
quantum, but not thermal, order by disorder. 

VII. SUMMARY AND OUTLOOK 

In this paper, we have considered quantum fluctua- 
tions of frustrated magnets in a magnetic field in a semi- 
classical framework, focussing on the issue of dynamical 
symmetries. In closing, we remark on some connections 
between this work and experiment, as well as other the- 
ories of frustrated magnets. 

A. Energy scales and experiment 

Even though all energy scales in the semiclassical treat- 
ment are set by JS, the prefactors which describe at 
what temperature different effects will become visible dif- 
fer greatly. Although none of these prefactors are large, 
they differ considerably in their smallness. 

The tendency towards collinear ordering (Fig. |2J - 
where possible - is the largest effect. With four bonds per 
unit cell having non-collinear spins, and three spins per 
unit cell, a zero-point energy per spin of 0.069 JS cos 2 9 
implies an effective coefficient of the biquadratic term 
of about JS/20. This might just be visible in the most 
frustrated magnets. If we take a jarosite with a spin as 
large as S — 5/2, there is a factor of 50 between JS 2 and 
JS/20; this is comparable to the suppression of the or- 
dering temperature due to frustration, which is typically 
also lowered by a factor of about 50 from the Curie- Weiss 
temperature^ 

Selection among those collinear states is weaker - the 
zero-point energy per spin differs by at most J5/50. To 
observe this effect, one requires not only the absence 
of other terms in the Hamiltonian (beyond the classi- 
cal nearest-neighbour exchange) on that energy scale. 
Rather, the system might not be able to equilibrate into 
the eventual ground state ensemble for reasons of 'topo- 
logical glassiness'fi^ that is to say it might not be able on 
experimental timescales to tunnel through the collinear 
anisotropy barriers. From the data plotted in Fig. [HJ it 
would nonetheless seem that the best place to look for 
such quantum order by disorder is near fields for which 
collinear configurations are available. Luckily, a breaking 
of a Z$ symmetry goes along with a clear signature in the 
specific heat, as has been observed elsewhere for kagome 
magnets with a Potts order parameter pi 

B. (Effective) gauge theories 

The general medial lattice construction A — > T makes 
the emergence of gauge structures in highly frustrated 
magnetism look rather natural - the spins for systems of 
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corner-sharing simplices quite naturally reside on links 
of the simplex lattice, A. A particularly striking mani- 
festation of the tendency of frustrated magnets to self- 
generate gauge descriptions is provided by dipolar spin 
iceM 

The classical ground-state constraint is a simple local 
constraint on the sites of A. For the theories considered 
here, it is a U(l) constraint on the lattice level, which as 
we have shown can be tuned by the application of a mag- 
netic field (the magnitude of the at each site of A is fixed 
by h). The effective long- wavelength gauge theories de- 
pend on further details of the lattice under consideration, 
such as bipartitenessii&M However, it is not so surprising 
from the current perspective that applying a field h by 
itself does not change the effective theory drastically, 35 
so that for example the different possible n-dimer models 
display very similar phase structures'^ 

In the context of the gauge theories describing the cor- 
relations in such dimer models, it is natural to think of 
the spins as (gauge-invariant) electric fluxes residing on 
the links of A. However, in the gauge-like theories of 
the Henley type considered here, the collinear spins have 
a second, complementary role as an Ising gauge (rather 
than gauge-invariant electric) field, the magnetic fluxes 
$ of which (Eq. I5.5fl determine the zero-point energies. 

It is this 'dual' role of the spins in the current con- 
text which leads to many of the unusual features of the 
problems under consideration. This is yet another exam- 
ple of (otherwise familiar and well-understood) effective 
theories in frustrated magnets exhibiting very unusual 
supplementary structures. 
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APPENDIX A: EQUATION OF MOTION IN A 
FIELD 

A spin precesses in its effective (exchange+Zecman) 
field: 

Si = i[Si,H] =-Y,SjXSi + hxSi 

j:i 

Summing this over i G v(a) gives 

L a = ^ S{ = ^ ^ Sj x Si ~\~ h X L a . 



Using Saf3 — L a + L/3 — 2S a p, this becomes 
L a = - y^^Lg + Lp) x S a p + hx L a . 

As ^2p :a S a p = L a , 

L a = + xS a p + hx (l a + ^) 

=> {j t -\h*)l a = ^2S aP x Lp . 

As the I are small deviations from the ground state, 
this expression can be linearised by replacing the spins 
{Sap} by their ground state values {S a °p}: 

d-hx)l a = J2S^xL . (Al) 



APPENDIX B: SPECTRUM FOR Q = AND 
(3 = ^x73 PLATEAU CONFIGURATIONS 

Iterating Eq. 13.31 once to eliminate the degrees of free- 
dom on one sublattice of triangles (corresponding to one 
sublattice of the bipartite honeycomb lattice A) gives 

( w -^) 2 c = EEw 7 + - (bi) 

In the following, we define the unit vectors e as e\ = 
(1,0), e 2 = (-1/2,73/2) and e 2 = (-1/2,-^3/2). We 
have chosen twice the nearest-neighbour distance of the 
kagome lattice as our unit of length. 



1. Q = 

The resulting spectrum is 



u) = ±— ± \Ji + 2(cosgi — cosg2 — cos 93) , (B2) 

where qt — q 1 e\ + q 2 ej, and the two choices of ± can be 
made independently. 

The zero-point energy per spin is then the average 
E = (M)/2. We find it to equal 3 ^ 

E = 0.5484JS' . (B3) 

Note that the simple form uu = ±-| ± is deceptive, as 
the square root can be larger than h/2, so that l^-l-ti^l = 
h ^ |«i| + \iv 2 \. 
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2. Q = v / 3x v / 3 

To compute the frequencies in this case requires solving 
the characteristic polynomial of a 3 x 3 matrix. Denoting 
u) = ±t| ± %/3 + A, this polynomial is given by 

-A 3 + /A + 5 = 0, (B4) 

Here, / = 9-2c_, g = 4+2c 3 -2c_, with q t = q 1 e\ + q 2 e 2 i 
and 

c_ = cos(qi - q 2 ) + cos(q 2 - q 3 ) + cos(q 3 - q x ) 



c 3 = cos(3qi) + cos(3<7 2 ) + cos(3g 3 ) . 

The frequencies are cumbersome to write down, but a 
closed form expression can be obtained by substituting 
the above into the formula for solving the cubic equation 
Eci.lBl 

Finally, we find a zero-point energy of 

E = 0.5643 JS . (B5) 
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